Thermal transistor: Heat flux switching and modulating 
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Thermal transistor is an efficient heat control device which can act as a heat switch as well as 
a heat modulator. In this paper, we study systematically one- dimensional and two-dimensional 
thermal transistors. In particular, we show how to improve significantly the efficiency of the one- 
dimensional thermal transistor. The study is also extended to the design of two-dimensional thermal 
transistor by coupling different anharmonic lattices such as the Frenkel-Kontorova and the Fermi- 
Pasta-Ulam lattices. Analogy between anharmonic lattices and single-walled carbon nanotube is 
drawn and possible experimental realization with multi- walled nanotube is suggested. 
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I. INTRODUCTION 

Much attention has been devoted in the past years to 
the study of heat conduction in low-dimensional systems 
Such study is not only important for understanding 
the fundamentals of statistical mechanics, but also for 
potential applications in heat control and management 
in nanoscale devices. 

Indeed, various models for thermal rectifiers/diodes 
that allow heat to flow easily in one direction have been 
proposed @, H, 0, i, 1, 0, 11, H- Moreover, based on 
the phenomenon of the negative differential thermal resis- 
tance (NDTR) observed in nonlinear lattices, a theoreti- 
cal model for "thermal transistor" is also proposed [l(|, 
which allows us to control heat current (due to phonons) 
by adjusting the temperature of the gate terminal (or 
called "gate temperature"). This is similar to the elec- 
tronic transistor which controls the electric current by 
adjusting the gate voltage. We should point out that our 
"thermal transistor" is different from the recent "heat 
transistor" , which controls the heat flux (caused by elec- 
trons) via adjusting the gate voltage Most recently, 
based on the thermal transistor, the different thermal 
logic gats have been also realized theoretically [12], which 
means that the phonons, the heat carrier, can also be 
used to carry information and processed accordingly. 

The purpose of the current paper are two folds. First, 
we shall discuss how to improve significantly the switch- 
ing efficiency of ID thermal transistor; Second, we shall 
extend the study from ID model to 2D which is closer to 
laboratory fabricated materials. 



II. ONE-DIMENSIONAL MODEL 

The ID thermal transistor model consists of three 
segments S, D, and G and each segment is a Frenkel- 
Kontorova (FK) lattice [HI, Qjj]. The model configura- 
tion, which is shown in Fig. 1(a), is similar to the one 
proposed in Ref.fioj except for an additional interface 



coupling ki n ts- In each segment, we regard the particle 
that is coupled to heat bath as the first particle and the 
interface particle as the last particle. The total Hamilto- 
nian of the ID model writes 



H = H S + H D +H G + H„ 



(1) 



and the Hamiltonian of each segment can be written as 
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with xw,i and Pw,i denote the displacement from equi- 
librium position and the conjugate momentum of the i th 
particle in segment W, where W stands for S, D or G. 
The parameters k and V are the harmonic spring con- 
stant and the on-site potential of the FK lattice. We cou- 
ple the last particle of segment S, D, and G to particle O 



via harmonic springs. Thus 
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and each segment 



Fixed boundaries are used, i.e., xw,o- 
consists of 10 particles. 

In our simulations, we use the Nose-Hoover heat baths 
[IH and we set the mass of each particle and the Boltz- 
mann constant to unity unless otherwise stated. In the 
nonequilibrium stationary state, the heat flux, which is 
the rate of energy transport, is constant along each seg- 
ment and the local equilibrium temperature is given by 
the time average of the square of particle's velocity per 
degree of freedom. 

We show in Fig. 1 that the model described by Hamil- 
tonian (1) can display different heat controlling features 
depends on the applied gate temperature, Tq. In par- 
ticular, the model can act as a heat switch, namely, the 
system can be switched from a heat insulator to a heat 
conductor by just adjusting the gate temperature. As 
demonstrated in Fig. 1(b), when Tq is at 0.01, 0.08, 
and 0.14, Js—Jd, thus Jq—0, namely, the control termi- 
nal does not provide any current to the system, the sys- 
tem, however is switched from an insulating (off) state, 
to a semi-conducting (semi-on) state, and conducting 
(on) state. At these three points, Jd,s are 3.8xl0~ 7 , 
1.8xl0~ 4 , and 3.2xl0 -4 , respectively. The ratio of Jd,s 
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FIG. 1: (a) Configuration of the thermal transistor. Pa- 
rameters of the system are (b) Thermal switch: Ts=0.01, 
Vs=5.00, fc s =0.35, fc ints =0.30; T D =0.20, Vb=1.00, fcz,=0.20, 
fci„ t£ .=0.05; V G =5.00, fc G =1.00, fc IJltG =0.20, and Vb=5.00. 
(c) Thermal modulator: all parameters remain the same as 
in (b) except with fci ntG =0.05. 



at the "on" state to that at the "off' state is considered 
as the switching efficiency which is about 800. 

We further demonstrate that the model can act as a 
heat modulator as well in which Jd.s can be continu- 
ously modulated by the gate temperature. To show this, 
we recall that for any two weakly coupled lattices, the 
heat current depends among others, on the strength of 
the interface coupling used [3j,ll6j]. Thus, to have Jd~Js, 
we can reduce the strength of ki nt o so that the heat 
flow through the control terminal decreases. Apparently, 
kintG should not be too small otherwise To is not well- 
controlled by Tq. The numerical result is shown in Fig. 

He). 



In our next configuration, we replace segment D with 
a FPU-/3 lattice whose phonon spectrum is contributed 
by low frequency vibrations (acoustic phonons) only. The 

Hamiltonian is given by H=^2iLi ^i^ kFFU '''' 
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4 (xi — Xi—i) , with being the anharmonicity. We show 
in Fig. 2(a) that such a setup reduces the heat cur- 
rent at "off" state greatly but affects the heat current 
at "on" state only slightly with the switching of about 
1300, where J D ,s equals to 2.15xl0 -7 and 2.82xl0 -4 at 
Tg=0.01 and 0.14, respectively. It can be noticed that for 
this and the former lattice choice in the present model, 
Jd.s at the "off" and "on" state is about 10 -7 and 10~ 4 , 
respectively, with a switching of about 10 3 , this is a sig- 
nificant improvement compared with previous one [lfj. 
The high switching is mainly a result of the significant 
reduction of Jd.s at the "off" state, which is basically 
due to the lower To at "off" state ( in the previous work 
it is designed at 0.04). Moreover, we found that using 
lattice with lower heat resistance as terminal D is ideal 
for achieving a large Jd- In Fig. 2(b), we show the heat 
modulating effect with the same setup. Since the demon- 
strated features of the transistor rely very much on an 
efficient control of the interface temperature, it would be 
interesting also to study how the property of the control 
terminal can affect its functions. In Fig. 3, we show that 
when we use the FPU-/? lattice as the control segment, 
due to the large heat resistance in the segment, temper- 
ature To does not well follow To thus heat currents in 
segment S and D do not sensitively depend on To in the 
low temperature regime. Only when Tq> 0.08, signif- 
icant increment of currents and continuous modulation 
of Jd,s are observed. In other words, Jd.s can be even 
more precisely adjusted. 

The results presented so far have shown that by choos- 
ing the lattice in different segment appropriately, one can 
improve the switching and modulating efficiency signif- 
icantly. This message is very important and useful for 
optimization and experimental realization of the thermal 
transistor. In the following, we shall extend the study of 
the thermal transistor to 2D. Such extension is important 
in application, as it takes into account the interaction be- 
tween longitudinal and transverse motions and therefore 
is closer to real material such as thin film etc. However, it 
is not our intention here to look for optimum switching/ 
modulating condition for this model. 



III. TWO-DIMENSIONAL MODEL 



The displacement from equilibrium position and the 
momentum of a particle in 2D lattice are labeled as 
q = {q x ,q v } and p = {p x ,p v }, respectively. The relative 
separation between two particles is Ar= | r \ = \ q{ — qi \ . The 
Hamiltonian of our model consists of the Hamiltonian of 
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FIG. 3: The thermal modulator in which the FK lattice 
is used as segment 5 and the FPU-/3 lattices are used as 
segment D and G. The system parameters are Ts=0.01, 
V s =5.00, fc s =0.35, fc lnts =0.60; T D =0.20, k D ^ FPU =0.20, 
/3 D =0.05, k in tn=0.05; k G -FPu=0-20, /3 G =2.00, k in tG=0.20, 
and V o =5.00. 



FIG. 2: (a) Thermal switch: a FPU-/3 lattice with 
fc.o(=fc.FPE/)=0.20 and /3=0.05 is used as segment D. Other 
parameters remain the same as in Fig. 1 (b). (b) Thermal 
modulator: all parameters remain the same as in 2 (a) except 
with £ws=0.25, £wd=0.055, and fc;„ t G=0.04. 
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and the Hamiltonian of the interface: 
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Uo{qO,j) + Vints(ArN s , x ,j;O t j) 



+Vi n tD(ArNo, x ,j;0,j) + Vi n tG(ArN GiXt j;0,j)]- (3) 

Similarly, all segments are coupled to each other via 
harmonic springs. In our model, segments S and G are 
FK lattices and segment D is FPU-/3 lattice. Thus, 
V s/G/mtw (Ar)=lAr\ V D (Ar)=|Ar 2 + f Ar 4 , (7 D =0, 

and Ug/Q/o(q) — -^yf- cos2nq x cos2irq y . The harmonic 
spring constant k and the on-site potential V may take 
different values in different segments. Periodic and fixed 



boundaries are used in the Y (transverse) and X (lon- 
gitudinal) direction respectively, namely, qi w ,o=qi Wt N Y 
and qo w ,j={Q, j) for ...,Ny}- Thus the topologies 

of the segments, as shown in Fig. 4(a), are cylinders. 
Nose-Hoover heat baths are coupled to the ends of seg- 
ment (i=l)- In the stationary state, temperature gra- 
dients are set up in the X direction and thus only heat 
currents through longitudinal links, Ji w {=}2^Zi Jiw -j) 
need to be taken into account, which are constant along 
each segment. 

In Fig. 4(b), we show the heat switching for the 2D 
thermal transistor model. It can be seen that when is 
at 0.01 and 0.07, Js=Jd and Jg—0- At these two points, 
Jd.s ar e 8.43xl0 -5 and 3.38xl0~ 4 respectively and thus 
the switching efficiency is about 4. This is a significant 
reduction when compared with the ID case and this sce- 
nario is also quite similar to the 2D thermal rectifier []|. 
We show in Fig. 4(c) that when T$ is increased, the 
difference between Jg and Jd is reduced, which gives a 
better modulation of Jd,s by the gate temperature. It is 
worth mentioning that in 2D case, the crossover energy, 
particles above which are able to overcome the on-site po- 
tential barrier is much lower than that in ID case. Thus 
the lattice parameters should be re-designed, otherwise 
the switching efficiency is reduced significantly as shown 
above. 
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FIG. 4: (a) The prototype of a two-dimensional thermal tran- 
sistor. Parameters used are (b) Thermal switch: Ts=0.01, 
&s=0.20, Vs=5.00, fc, ltS =0.50; Zb=0.17, k D -Fpu=1.00, 
/3=3.00, fc mto =0.035; k G -FPu=l-00, Vg=6.00, k int G=2.50, 
and Vb=5.00. (c) Thermal modulator with T s =0.015. Here, 
N D ,x=N G ,x = W, N s ,x=9, and ATy=4. 



IV. DISCUSSION AND CONCLUSION 



merical simulations have demonstrated that such a pro- 
totype can work as a heat switch and a heat modulator 
- two basic functions of a thermal transistor. 

In our numerical simulations, we have used dimension- 
less units for the sake of computations. However, all di- 
mensionless units are related to the dimensional physical 
variables through dimensional scaling and consequently, 
the physical (real) temperature T r is related to the nu- 

muj 2 2 

merical temperature T through T r = — ]^~T where ks is 
the Boltzmann constant (please refer to [13] for details). 
Suppose we use single- walled carbon nanotube as our ex- 
perimental prototype for the thermal transistor model, 
namely with the mass of particle m, lattice constant a, 
and natural oscillation frequency tu correspond to the 
real values of such a material P_9(, T=0.01 in our simu- 
lations corresponds to physical temperature T r of about 
100 K. Similar analysis for heat flux gives J r (^muj^a 2 J) 
~10~ 5 J. Therefore, with J ~10~ 4 , J r is about 10" 9 W 
(J/s) which is smaller than that one in Ref. 20]. The 
smaller heat flux obtained in our model is not surprising. 
It is mainly due to the interacting with external on-site 
potential 2l| and the interface thermal resistance arises 
in coupling dissimilar nonlinear lattices Q , which is sim- 
ilar to the case of electronic diode where the electrical 
current are greatly reduced by the P-N junction. For 
a homogeneous FPU lattice with similar parameters, J r 
can be about 2 orders of magnitude higher. 

Even though a real working prototype for solid state 
thermal transistor is still an open issue, we believe that 
the above model or its variants will be realized sooner or 
later. Moreover, a solid state thermal rectifier using nan- 
otubes has been experimentally demonstrated [22| (see 
also Ref[23| for other variants). More recently, a tun- 
able and reversible thermal link using multi-walled car- 
bon nanotubes (MWCNTs) [2J] were successfully demon- 
strated whereby the thermal conductance can be con- 
trolled mechanically by displacing the outer shell with 
respect to the inner one. This is a step forward to the re- 
alization of thermal transistor. In fact, in the heat control 
models, be it a thermal rectifier or a thermal transistor, 
the key factors are: (1) the broken of the spatial symme- 
try of terminal, namely they are two different materials; 
(2) nonlinearity in each segment. Therefore, the possi- 
ble material to be used for the thermal transistor might 
be multiwalled nanotubes. As the recent experiment [241 ] 
shows that the outer wall can be used to control the ther- 
mal conductance of the nanotube, thus it can be used as 
the controlling terminal G. 



In summary, we have provided several ways in im- 
proving the switching efficiency in ID thermal transis- 
tor, namely, one can change different lattice in differ- 
ent segments, and/or by adjusting the coupling strength 
at the interface. Our study has been extended to two- 
dimensional lattices and a prototype of 2D thermal tran- 
sistor - a tube-like model have been proposed. Our nu- 
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